Looking for a light Higgs boson in the overlooked channel 
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The final state obtained when a Higgs boson decays to a photon and a Z boson has been mostly 
overlooked in current searches for a light Higgs boson. However, when the Z boson decays leptoni- 
cally, all final state particles in this channel can be measured, allowing for accurate reconstructions of 
the Higgs mass and angular correlations. We determine the sensitivity of the Large Hadron Collider 
(LHC) running at center of masses energies of 8 and 14 TeV to Standard Model (SM ) Higgs bosons 
with masses in the 120 — 130 GeV range. For the 8 TeV LHC, sensitivity to several times the the 
SM cross section times branching ratio may be obtained with 20 inverse femtobarns of integrated 
luminosity, while for the 14 TeV LHC, the SM rate is probed with about 100 inverse femtobarns of 
integrated luminosity. 



Introduction - The search for the Higgs boson is en- 
tering a critical phase. Data collected at the LHC rules 
out the SM Higgs boson for a wide range of masses 
and may suggest a Higgs boson with mass near 125 
GeV [TJ [5]. Searches for a light SM Higgs in the still- 
relevant mass window rely primarily on the 77 channel, 
though the WW* — > 2l2v channel and the golden chan- 
nel, ZZ* — > 4£, are also important. 

So far very little attention has been given to the 
Z7 — > ££j channel [3], although its event rate is compa- 
rable to that of the golden channel for a light SM Higgs 
boson. Nevertheless, this channel has the advantage that 
all final state particles can be measured well, which car- 
ries several important implications: 1) the Higgs mass 
could be measured from the total invariant mass spec- 
trum, 2) the spin of a putative signal can be determined 
by studying angular correlations [4 , and 3) the sepa- 
ration of signal from background can be facilitated by 
employing full kinematic information, potentially allow- 
ing searches with enhanced sensitivities. For the golden 
channel in ZZ* — > A£ the above questions have been stud- 
ied extensively [5HZ] > but we are not aware of any detailed 
studies for the Zj channel. 

Measurements of all four Higgs decay modes into elec- 
troweak bosons are in fact very important in determining 
the electroweak quantum numbers of a putative Higgs 
signal [8j. Furthermore, an electroweak singlet scalar 
could easily have a branching fraction in the Z"f mode 
that is orders of magnitude larger than the SM expecta- 
tion [5] which provides an important additional incentive 
for studying this channel. 

Kinematics - The kinematics of Zj — > ££j is de- 
scribed by three angles, 0, and (f>, where O may be 
taken to be the angle describing the production of the 
Z boson in the center of mass frame, and 9 and (f> are 
the angles that describe the decay of the Z to leptons, as 
defined in more detail in 0. To accomodate events with 
jet radiation in the final state, we use the momentum of 
the Zj system in the lab frame, rather than the beam 
axis, to define 9. 

The dominant irreducible background to the Higgs sig- 
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FIG. 1: Feynman diagrams contributing to qq 
shown in (a) and (b). 



li'y are 



nal arises from initial state radiation (ISR) and final state 
radiation (FSR) from Drell-Yan production of a Z boson; 
the diagram describing this process is shown in Fig. [I] (a) 
and (b). The invariant mass of the Zj system from FSR 
events is close to the Z boson mass, so this background 
is removed efficiently by imposing m« 7 > 100 GeV, and 
we can focus on the ISR diagram and the corresponding 
it-channel diagram for the rest of this analysis. 

The signal and background cross sections were com- 
puted using the helicity basis in [TO]. We now discuss 
some qualitative features of these differential cross sec- 
tions, in particular the dependence of the signal and 
background processes. 

In the signal case, angular distributions follow from the 
fact that the Higgs is a scalar particle, and hence only 
the decay angle 9 has a nontrivial distibution: 



1 



da 



N d cos 9 d cos 9 1 



= (1 + cos 2 



(1) 



For the background distributions, the non-vanishing he- 
licity combinations arc (Ai,A2) = (±, =p), (0, ±), and 
(±,±). The production angular distribution exhibits a 
collinear singularity at cosG = ±1, which is seen by ex- 
amining the t-channel propagator in Fig. [T] (a) , 

1 = 1 (2) 

(kq- Pl ) 2 2S 9 -S 7 (l-cos6) ' v ; 

while the u-channel propagator gives the collinear sin- 
gularity at cos = —1. Thus the production angu- 
lar distribution for the background process is peaked at 
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FIG. 2: Signal (red, solid) and background (blue, dashed) distributions in cos©, </> and cosd, with Vs = my, = 125 GeV. 



cos O = ±1, producing forward and backward photons. 
The singularity is removed by the px cuts on the photon 
and leptons. Explicit calculations lead to 
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with 
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where gz,(i) and 9rm are the Z couplings to left- and 
right-handed quarks (leptons). 

In Fig. [2] we show the distributions in cos©, <fi, and 
cos 9 for a 125 GeV Higgs boson and a background pro- 
cess dd — > Z^j at Vs = 125 GeV at the parton level. 
These are modified after including the effects of parton 
distribution functions (PDF) and detector acceptance 
and isolation cuts. In particular, we note that cos is 
directly connected to the photon px through 
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The cos distribution in Fig.[2]therefore implies that the 
p 7 T distribution is peaked at zero for the background and 
(to 2 — m 2 z ) / (2rrih) for the signal. However it also follows 
that once a cut on p 7 x is imposed, very little additional 
sensitivity can be gained from the cos distribution. 

Analysis and Results - We perform Monte Carlo 
simulations to obtain projections for the sensitivity of 
this channel at the LHC using various analyses. We con- 
sider Higgs masses of 120, 125, and 130 GeV. Our simula- 
tions are specific to the 8 and 14 TeV LHC. The existing 
7 TeV data has very little sensitivity in this channel such 
that we do not report those results here. 



To perform these Monte Carlo studies, we generate at 
least 50,000 events for each signal and background pro- 
cess using MadGraph 5 [TT]. The Higgs coupling to 
gluons and the hZj vertex are implemented as effective 
dimension five operators using the HEFT model provided 
by MadGraph 5 and the FeynRules package. For 
both signal and background, the processes pp — > and 
pp — > Z'y + lj are generated, using the MLM matching 
scheme |13j implemented in MadGraph 5 and interfaced 
with Pythia 6 [Hj, with a matching scale of 25 GeV. 
Events are then passed to PGS 4 [TS] using the CMS pa- 
rameter card, to model detector acceptance and smearing 
effects. 

Since the energy and momentum resolution is crucial 
for this analysis, we have compared the invariant mass 
resolution obtained from PGS 4 with the one that is ob- 
tained when smearing parton level events by hand using 
the CMS detector parameters [TB], and found that they 
agree in general. 

We demand that each lepton or photon has 



\t}\ < 2.5 and p T > 15 GeV. 



(7) 



The smearing results in the broadening of the lineshape 
in the total invariant mass of the Z"f system, to^^ 7 , for 
the signal events. Therefore, before performing more de- 
tailed analyses, we perform an invariant mass cut; de- 
manding that the invariant mass of the Zj system be 
within 5 GeV of the mean invariant mass of the Z7 sys- 
tem in signal events. It is worth emphasizing that since 
subsequent analyses will effectively reduce the range of 
invariant mass considered, the specific details of this ini- 
tial cut does not have a strong effect on the final value 
of S/\ /f B obtained. Note that this cut also effectively re- 
moves the background coming from FSR radiation that 
is characterized by to« 7 ~ Mz- 

To determine the expected number of signal events at 
the 14 TeV LHC, we obtain the inclusive Higgs produc- 
tion cross section from [T7]. For the 8 TeV LHC, we 
use the values given in [IB] . The branching fraction for 
h — > Zj is found using HDECAY [19], while we use the 
PDG value (6.73%) for the branching fraction for a Z 
decaying to leptons [5D]. The background cross section 
is found by using MCFM g3[52! with FSR photon radi- 
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8 TeV LHC, 20 Inverse Femtobarns 
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FIG. 3: Exclusion limits at the 95% confidence level on the 
Higgs production rate times branching fraction to Z'y at the 
8 TeV LHC with an integrated luminosity of 20 ficT 1 . The 
green (yellow) band is the 1(2) a contour. The solid red line 
corresponds to the SM expectation. 



Higgs Mass 


Signal (fb) 


Backg. (fb) 


S/VB (20 fb" 1 ) 


120 GeV 


0.38 (0.45) 


32. (110) 


0.30 (0.19) 


125 GeV 


0.61 (0.74) 


30. (100) 


0.50 (0.33) 


130 GeV 


0.66 (0.86) 


23. (89.) 


0.62 (0.41) 



TABLE I: The signal and background cross sections, as well 
as the significance after an optimal cut on the discriminant 
in Eq. in the invariant mass only analysis at the 8 TeV 
LHC. In the parenthesis we also show the corresponding values 
for all events passing the pt and geometric acceptance cuts 
and which are within an invariant mass window of 10 GeV 
centered on the Higgs mass, as described in the text. 



Higgs Mass 


Signal (fb) 


Backg. (fb) 


S/y/B (100 fb" 1 ) 


120 GeV 


0.83 (1.0) 


36. (180) 


1.2 (0.78) 


125 GeV 


1.3 (1.6) 


37. (160) 


2.0 (1.3) 


130 GeV 


1.7 (2.1) 


40. (140) 


2.7 (1.8) 



TABLE II: Same as Tab. for the 14 TeV LHC, with a 
luminosity of 100 fb . 



ation turned off. 

We perform three analyses, two of which are multi- 
variate. The multivariate discriminants we use are based 
on the matrix element of the signal and background pro- 
cesses. In the context of a maximum likelihood analysis 
such a discriminant was used in the discovery of the single 
top production in [53] . For simplicity we use a cut-based 
approach to determining our sensitivity using these mul- 
tivariate discriminants. 

We construct a discriminant using the fully differential 
cross sections computed for the signal and background 
processes to quantify the relative probability of a partic- 
ular event being signal-like or background-like. We then 
determine an optimal cut on the discriminant to maxi- 
mize the value for S/\/~B. In one analysis, we include 
PDF weights for the leading initial state for signal or 
background events (gg or qq respectively) . In the second 
multivariate analysis, we do not include a weight from 
PDFs. Labelling the signal and background differential 
cross sections by s(f2) and b(Sl), respectively, we consider 
the quantity 

Here, ft — {x\, x<i, s, rn e g, O, 9, 4>} is the complete set of 
kinematic observables characterizing each event. When 
evaluating D on a sample of pure signal events the dis- 
tribution is peaked toward 1 while it is peaked toward 
for a pure background sample. For each Higgs mass, 
a cut on D is determined by maximizing S/ \B of the 
events passing the cut. One advantage of using the mul- 
tivariate discriminant in a cut-based approach is that the 
relative normalization of the signal and background cross 
sections does not affect the final significance computed 
using S/yB. The drawback, on the other hand, is that 
we lose those events not passing the cut, which would not 



be the case if the signal and background matrix elements 
were used to construct the likelihood directly. 

Our multivariate discriminants use the parton-level dif- 
ferential cross section except for the Higgs propagator, as 
for the Higgs masses considered the Higgs width is much 
narrower than the experimental resolution. In princi- 
ple, one can deal with this issue by using transfer func- 
tions for the lepton momenta. We take the simpler ap- 
proach of weighting each event with a Gaussian invariant 
mass distribution that is centered at the average invariant 
mass for signal events. The width used in this Gaussian 
weighting is found by scanning (in 20 MeV increments) 
over potential values, from 100 MeV to 5 GeV, and select- 
ing the value which maximizes the sensitivity of the anal- 
ysis. The third analysis uses the same Gaussian invariant 
mass weight, but no other kinematic information about 
the events. While one would expect a loss of sensitivity, 
this approach has the advantage of being less sensitive to 
higher order corrections that could modify the angular 
distributions that enter the multivariate analyses. 

We find the best values for S/ \[B from the analysis in 
which the full differential cross sections and PDF weights 
are used. However the sensitivity from this analysis is 
only ~ 1% larger than that obtained from the invariant 
mass only analysis. The smallness of this increase in sen- 
sitivity is due to the fact that the relatively hard p 7 T 
cut leaves us without much additional sensitivity to O, 
and the other angular variables are not as sensitive, espe- 
cially given geometric acceptance and finite momentum 
resolution. We therefore quote results using the invari- 
ant mass only analysis, as they should be more robust 
with respect to systematic uncertainties. In particular, 
the m« 7 distribution is unaffected by jet radiation, so 
that corrections to the jet multiplicity and momentum 
distribution, which is only simulated to leading order in 
our analysis, will not reduce the sensitivity. 

The signal and background cross sections after the op- 
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14 TeV LHC, 100 Inverse Femtobarns 
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FIG. 4: Exclusion limits at the 95% confidence level on the 
Higgs production rate times branching fraction to Z'y at the 
14 TeV LHC with an integrated luminosity of 100 flT 1 . The 
green (yellow) band is the 1(2) a contour. The solid red line 
corresponds to the SM expectation. 

timal cut on D from this invariant mass only analysis are 
listed in Table [I] for various Higgs masses at the 8 TeV 
LHC. The expected significance with 20 fb -1 integrated 
luminosity is also provided. Tabic [TT] shows analogous in- 
formation for the 14 TeV; here the expected significance 
with 100 fb _1 is shown. 

In the absence of any signal, we have also considered 
the expected exclusion limit on the Higgs production rate 
in the gluon fusion channel using the CL S method [21] 



with 20 fb 1 of integrated luminosity for the 8 TeV LHC 
in Fig.[3]and for the 14 TeV LHC with 100 ftr 1 in Fig. [4] 

Conclusions - We have considered the possibility of 
searching for a light Higgs boson in its decays to ££7 final 
states via Z7. This branching ratio is known precisely in 
the SM, and deviations from this rate are unambiguous 
signals of new physics that couples to the Higgs boson, 
or could even signal the presence of a Higgs imposter [5] . 

We have performed a detailed Monte Carlo study for 
the 8 and 14 TeV LHC. We find that branching ratios 
for the Higgs decay to Z7 of several times the SM rate 
are probed at 8 TeV with 20 fb -1 , while the SM rate is 
probed at the 14 TeV LHC with 100 fb -1 . For Higgs 
masses of 125 GeV and above, a measurement of the 
Higgs branching ratio to Z7 is in reach of the 14 TeV 
LHC. We hope this work inspires experimental efforts in 
this particular search channel. 
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